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, ^ ■ Abstract 
^ . 

We examine the spectrum of super symmetric particles predicted by grand 
unified theoretical (GUT) models where the electroweak symmetry breaking is 
accomplished radiatively. We evolve the soft supersymmetry breaking parame- 
^ ■ ters according to the renormalization group equations (RGE). The minimization 

I of the Higgs potential is conveniently described by means of tadpole diagrams. 

■ We present complete one-loop expressions for these minimization conditions, 

including contributions from the matter and the gauge sectors. We concentrate 
on the low tan (3 fixed point region (that provides a natural explanation of a 



I large top quark mass) for which we find solutions to the RGE satisfying both 

experimental bounds and fine-tuning criteria. We also find that the constraint 



from the consideration of the lightest supersymmetric particle as the dark mat- 
ter of the universe is accommodated in much of parameter space where the 
lightest neutralino is predominantly gaugino. The supersymmetric mass spec- 
. trum displays correlations that are model-independent over much of the GUT 

parameter space. 

X 

a • 1 Introduction 



Why should one be interested in supersymmetry? Until recently, the reasons have 
been principally theoretical. Supersymmetry (SUSY) is a beautiful extension of the 
Poincare symmetry with new dimensions of space and time that explain the existence 
of fermions||lI]. It solves the hierarchy problem of widely separated electroweak and 



grand unified scales through cancellations among diagrams that give quadratically 



divergent Higgs boson mass corrections. Moreover supersymmetry may be a necessary 
consequence of string theory. 

The recent upswing in interest in supersymmetry derives from high precision mea- 
surements of Standard Model (SM) parameters at LEP. Renormahzation group evolu- 
tion with minimal SM particle content of the SU(3), SU(2), and U(l) couplings from 
converge at a single high scale, in contradiction with the prediction 
of the SU(5) grand unified theory (GUT). However, with the minimal particle content 
of supersymmetry included, the evolution is in excellent agreement with LEP data 
and suggests a grand unified scale at Mq ^ 2 x 10^^ GeV and effective SUSY mass 
scale within the range < M^jjgy < 1 TeVp|. Encouraged by this success, the evo- 
lution of Yukawa couplings is also being vigorously pursued, with Yukawa unification 
constraints such as Xb = Xt at the GUT scale 0. While the unification of gauge and 
Yukawa couplings is an extremely attractive feature, the existence of supersymmetry 
will only be confirmed when new particle states are seen directly and the associated 
R-parity conservation or violation is tested in the production and decays of these 
supersymmetric particles. 

The idea of a radiative breaking of the electroweak symmetry is an old but still 
popular one[|]-[|l3l. It is very attractive to explain the breaking of the electroweak 
symmetry through large logarithms between the Planck scale and the weak scale 
For the radiative corrections to be strong enough to drive a Higgs boson mass-squared 
parameter negative (thus breaking the electroweak symmetry), a Yukawa coupling of 
that Higgs boson must be large at the GUT scale. With the top quark mass large 
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{nit > 100 GeV), the SUSY GUT unification can naturally explain the origin of the 
electroweak scale. A heavy top is required to drive one of the soft-supersymmetry 
breaking parameters (a Higgs doublet mass) negative. Today we know the top quark 
mass is large and that the top has a large Yukawa coupling. There is a relation- 
ship between the electroweak scale and the top quark Yukawa coupling through the 
RGEs; consequently the radiative symmetry breaking mechanism has important con- 
sequences for the supersymmetric particle spectrum. Indeed a large top Yukawa 



coupling is the motivation for the fixed point solutions 11^ advocated recently in the 



context of GUT theories[jT5|-|T8[. These solutions predict a linear relationship be- 
tween rrit and sin f3, given further constraints on the SUSY particle spectrum. 

There are at least two other motivations for supersymmetry. In the context of 
SUSY GUTs, the grand unification scale is raised sufficiently high to suppress proton 
decay to experimentally acceptable levels, when an additional R-parity symmetry is 
invoked. R-parity symmetry has an important consequence, providing the second ad- 
ditional motivation for supersymmetry - it implies that the lightest supersymmetric 
particle (LSP) is stable. It is now generally believed that baryonic matter is insuf- 
ficient to make up the total gravitationally interacting matter of the universe. The 
LSP provides a natural candidate for the (cold) dark matter of the universe, since 
the LSP is forbidden to decay into baryons by R-parity conservation. 

The greatest dilemma for supersymmetry is the unknown mechanism for super- 
symmetry breaking. Since no supersymmetric partners of the known particles have 
been observed, one must not only explain why supersymmetry is broken, but also why 
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it is broken in such a way that (almost) all of the supersymmetric partners are heav- 
ier than the known particles. Typically one skirts the issue of the exact mechanism 
of supersymmetry breaking, and parametrizes one's ignorance by introducing soft- 
supersymmetry breaking parameters which characterize the scale of the masses of the 
supersymmetric particles. In the context of certain models (minimal supergravity, 
no-scale supergravity, superstring-inspired, etc.) there are relations between these 
parameters at the GUT scale that appear as universal boundary conditions. The 
values of the parameters would be determined in some still unknown fundamental 
theory, but in practice the models are described as depending on a few undetermined 
soft-supersymmetry breaking masses which, when evolved in the well-known effective 
theory below the unification scale, generate consequences for the masses and couplings 
of the SUSY particles that can be compared to experiment. 

The supersymmetry breaking in such models occurs through a hidden sector at 
some scale intermediate between the electroweak and Planck scales. For example, 
there may exist an asymptotically free gauge theory in the hidden sector in which 
gaugino condensation or some other nonperturbative mechanism occurs at a high 
scale HgjjgY related by a geometric hierarchy to the Planck and electroweak scales 
{Mg^gY ~ H'sijgY / ^p^^)-i thus breaking supersymmetry in the observable sector near 
the electroweak scale Mg^gy ~ M^^. Since this breaking of supersymmetry occurs 
in the hidden sector which is coupled only gravitationally to the observable sector, 
the supersymmetry breaking is suppressed relative to the Planck scale by factors of 
Hg^jgy/Mp. This mechanism gives rise to an effective theory in the observable sector 
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with softly broken global supersymmetry, giving the scale of the soft-supersymmetry 
breaking parameters near the electroweak scale. How the effective theory parameters 
are determined and related to each other depends on the details of the supersymmetry 
breaking scenario. 

2 Soft-supersymmetry breaking parameters 

Retaining only the dominant Yukawa couplings Aj, A^ and A^, the superpotential* is 
given in terms of the superfields by 

W = (XtQ'Hif + XbQ'Hib' + XrL'Hir' + fiHlHi) , (1) 

where Q = {t,b), L = (r, Ut-) and Hi = (ifj*, iff) and Hi = {H2, H2) and Cjj with i, 
J = 1, 2 is the antisymmetric tensor in two dimensions. The Yukawa couplings are 
defined by 

At — — : — ^ , Ab — ^ , — o ' l^J 

vsmp V cos p V cos p 

where tan /3 = is the ratio of the vacuum expectation values of H2 and H^. The 
II term in the superpotential contributes to the Higgs potential which at tree level is 

Vo = {ml^+ii')\Hi\'' + {ml^ + ii^)\H2\' + rnl{eijHi'H2^ + h.c.) 
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+ -{9' + 9") \Hi\'-\H2\' +-9'\Hl*H^\\ (3) 
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where mff_^, rriff^, and ms are soft-supersymmetry breaking parameters. We shall 
define as usual the soft Higgs mass parameters 

mj = + /i^ , (4a) 



* Caution: Our convention for the sign of ^ (and also the sign of the trihnear couphngs A intro- 
duced below) differs from some authors. 



ml = mjj^ + /i^ . (4b) 

Of the eight degrees of freedom in the two Higgs doublets, three (G^, G°) are absorbed 
to give the W"^ and Z masses, leaving five physical Higgs bosons: the charged Higgs 
bosons H^, the CP-even Higgs bosons h and H, and the CP-odd Higgs boson A. 
There are soft-supersymmetry breaking gaugino mass terms 

^M^BB + ^M2W^W'' + ^MsW , (5) 

for the bino B, the winos (a = 1, 2, 3), and the gluinos {h = 1,...,8). 
Corresponding to each superpotential coupling there is a soft-supersymmetry breaking 
trilinear coupling 



[XtAtQ'Hif + \,A,Q'Hih' + KArVHif' + fiBHim) (6) 
and soft squark and slepton mass terms 

+ Ml[flfL + ylyL] + Mlf},fR. (7) 
The RGE for the soft-supersymmetry breaking parameters are given in the appendix. 



and the RGE for the gauge and Yukawa couplings are summarized in Ref. JTS . 

An interesting aspect of the supergravity breaking mechanism is the origin of the 
3 — 2 — 1 supersymmetry at low scales. Why is it the electroweak gauge group the one 
that is broken, and not QCD? Consider the renormalization group equations from the 
appendix for the scalar states H2, iji, and Ql retaining only the QCD gauge coupling 
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and the top Yukawa coupling At terms [[19 

16 
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where = M^^ + M^^ + Ml^^ + and t = In Q/Mg- The term is the means by 
which the mass-squares are driven to lower values as the scale decreases. Because the 
Higgs field is uncolored, the group theory factors allow to be driven negative with 
M^^ and remaining positive, thus breaking only the electroweak gauge group. 

According to conventional wisdom the squarks and sleptons have a universal soft- 
supersymmetry breaking mass at the unification scale. Then any deviations from 
degeneracy at the SUSY scale are suppressed by the associated quark or lepton mass, 
which is small except for the top squarks. The fiavor changing neutral currents 
(FCNCs) are thereby suppressed to an acceptable level. The universal boundary 
condition applies in minimal supergravity models with the canonical kinetic energy. 
Recently there has been some interest in relaxing this condition [pO|- p2| . 

Analytical expressions can be obtained for the squark and slepton mass param- 
eters when the corresponding Yukawa couplings are negligible (i.e. for the first two 
generations). For a universal scalar mass and gaugino mass mi at the GUT scale 

2 

(this condition need not apply in general in string theories), one has the relation 

3 

m^- = + ^ f(m\ + {T^j - ej sin^ ^^)M| cos 2/? , (9) 
1=1 

for the squark and slepton masses where the ft are (positive) constants that depend 
on the evolution of the gauge couplings 



fi 



1 - ^^ht 



(10) 



Here j is the SU(2) quantum number and ej is the electromagnetic charge of the 
sfermion. The bi are given in the appendix and Ci{f) is ^^j^ (0) for fundamental 
(singlet) representations of SU{N) and j^Y'^ for C/(l)y. The squark mass spectrum 
of the third generation is more complicated for two reasons: (1) the effects of the third 
generation Yukawa couplings need not be negligible, and (2) there can be substantial 
mixing between the left and right top squark fields (and left and right bottom squark 
fields for large tan /3) so that they are not the mass eigenstates. 

The gaugino evolution is particularly simple by virtue of their simple renormal- 
ization group equations; at one- loop order the gaugino masses parameters Mi, M2, 
and M3 scale in exactly the same proportions as do the gauge couplings so that 



Figure 1 shows a typical evolution of the soft-supersymmetry breaking param- 
eters. The characteristic behavior exhibited by the mass parameters are typical of 
renormalization group equation evolution. The colored particles are generally driven 
heavier at low Q by the large strong gauge coupling. The Higgs mass parameter m| 
is usually driven negative (at least for tan/3 not too small), giving the electroweak 
symmetry breaking. Assumed universal boundary conditions at the GUT scale yields 
correlations between the masses in the supersymmetric spectrum. 



rrig = Ms{mt) 



a2{mt) 



M2{mt) 



ai{mt) 



Mi{mt) . 



(11) 
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Evolution of sparticle masses 
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Fig. 1. An example of the running of the soft-supersymmetry breaking parameters 
for asiMz) = 0.120, mtimt) = 150 GeV, tan/3 = 10, mi = 250 GeV, mo = 100 

2 

GeV, and = 0, where the superscript G denotes the GUT scale. 



Fixed-point solutions to the RGE predict that the scale of the top-quark mass 
is naturally large in SUSY-GUT models but depends on tan/3. The prediction is 



that 15 



pole 

ml 



(200 GeV) sin 13 . 



(12) 



Note that the propagator-pole mass is related to this running mass nit{nit) 



77iP°'^ = mt{mt) 



1 + ^"3("^t) 
OTT 



(13) 
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3 One-Loop Contributions: Tadpole Method 

Although the tree level Higgs potential is not reliable for the purpose of analyzing 
radiative breaking of the electroweak symmetry]^], it provides a convenient starting 
point for our discussion. Recall the tree level potential Eq. (^. In supergravity 
models ml is related to B and by 

ml = Bfi . (14) 

When the neutral components of the Higgs doublets receive vacuum expectation 



values Vi and V2, the potential develops tadpoles. Inserting] 24] 



H,= (7=.(^^+:}+'^^^] (15a) 



Hi 



^2=( 1 ) (15b) 



into Eq. (§) one can identify 

Vtadpole = tiiJi+t2i^2, (16) 

where ti and t2 are (tree-level) tadpoles: 

^1 = {ml^+fi^)vi + BfiV2 + ^{g^ + g'^)viivl-vl) , (17a) 
t2 = {ml^+fi')v2 + BfiVi-^{g' + g'^)v2{vl-vl) . (17b) 

The minimum of the Higgs potential is determined by setting the first derivatives of 
the fields to zero, 

dVo dVtadpole 
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Therefore the tadpoles ti and ^2 must vanish at the minimum. With our normahzation 
of ipi and 4'2 (i-e. including the factor of in Eq. ( p.5a| ) and Eq. ( |15b| )), the W and 
Z masses are 



Ml 



-'- 2/ 2 I 2\ 

-^9 {v^+v^) , 



(19a) 
(19b) 



which implies f ^ + f | = f ^ = (246GeV)^. A particularly useful form of the mini- 
mization conditions is obtained by forming the linear combinations Ti and T2 of the 
tadpoles given by 

( T,\_ ( cos/5 -sin/5 

[t, )-[ sin/5 cos/5 )[h ) ' ^""^ 
where cos/5 = Vi/v and sin/5 = V2/V. From Eqs. ( |17aD and ( |17b| ) we have 

1 



To 



1 

V 
V 

1 r 

V 
V 



'm]j^ + /i^) cos^ /5 - {m% + /i^) sin^ (3 +\{g'^ + g''^)v^ cos 2/5 



+ "^^2 + 2/i^) sin 2/5 + 



(21a) 



(21b) 



We see that the rotation (^) through the angle /5 conveniently places all of the 
dependence on gauge couplings (D-terms) in Ti. Setting Ti = and dividing by 
V cos 2/5 yields the familiar tree- level condition 



1 



Ml 



m'jj_^ — m'jj^ tan^ /5 



2"'^ tan2/5-l 
Setting T2 = and dividing by v the other tree-level condition 

-Bfi = ^(m2,^+m|,^+ 2/^2) sin 2/5, 



(22) 



(23) 
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is obtained. Notice that the signs of B and /i are not determined by the minimization 
conditions (only the relative sign is known), giving rise to two distinct cases. 

We can extend the above technique to include one-loop contributions to the Higgs 
potential, deriving equations analogous to (P^ and ( P^D by setting to zero linear 
combinations of tadpoles rotated through the angle /3. The one-loop effective potential 
is given by 



Vo + 



(24) 



where Vq is the tree-level Higgs potential and 



1 



647r2 



rStr 



In 



(25) 



is the one-loop contribution given in the dimensional reduction {DR) renormalization 
scheme p5||. The supertrace is defined as Str/(A^^) = I]i Cj(— l)^'**(2sj + l)/(mf) 



where Ci is the color degrees of freedom and Si is the spin of the i*^ particle. To 
determine the minimum one must set the first derivatives of the effective potential to 



zero 



^ = ^ + ^Str 



M In ■ 



dip 



(26) 



We note that f{m'^) usually involves the mass eigenstates of the theory; one therefore 
ought to use the coupling of the Higgs fields to the mass eigenstates in tadpole calcu- 
lations in order to facilitate comparisons between minimization techniques. Evaluated 
at the minimum of Vi, tadpole contributions involve the coupling dM."^ /dip and the 
usual integration factor (in ^ — l) ; setting tadpole contributions to zero is 
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therefore equivalent to minimizing the potential. More generally, the nth derivatives 
of the effective potential are related to the diagrams (at zero external momentum) 
with n external lines; the minimization conditions at one-loop are obtained by calcu- 
lating diagrams with only one external line - the tadpoles PB|, P7| . 

In order to maintain the linear combinations in (|20|) for the tree level relations, we 
calculate with appropriate combinations of Higgs fields in the external Higgs line in 
the tadpole diagrams. The Feynman rules usually express these external Higgs lines 
as the physical Higgs bosons H or which are obtained from the Higgs fields z/'i, 11)2 
by a rotation by an angle a (in the opposite direction to the rotation (3 performed 
above) 

H \ _ ( cos a sin a \ / ipi 
hi \ — sin a cos a 1 \ ^2 

As with the tree-level tadpoles, we need to rotate the one-loop contributions by 

the same angle /3 in order to express the minimization conditions most simply. We 

therefore define the desired linear combinations J', J7j_ of Higgs fields 

( J \ _( cos f3 -sm(3 \ f ipi \ _ f cos(/? + a) - sm{/3 + a) \ f H \ 
\ J± J ~ \ sinP cosp ) \ip2 J ~ \ sin(/3 + a) cos{p + a) J \ h J ' 

(28) 

To include the one-loop corrections, we calculate the tadpole diagrams in Figure 2, 
and add the suitably regularized result to the tree-level results. 
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Fig. 2. The one-loop tadpole diagram. The loop consists of matter and gauge-Higgs 
contributions. 

This tadpole technique is not new, and is equivalent to procedures followed pre- 
viously. However it provides an alternate way of organizing the calculation and of 
understanding why the contributions have their particular form. Moreover, the an- 
alytical expressions obtained with the tadpole technique are often very useful, par- 
ticularly in certain regions of parameter space that are difficult to explore by simply 
minimizing the potential numerically (e.g. the low-tan /3 fixed-point region). 

The method of determining the minimization conditions at one-loop by calculating 
tadpoles is especially convenient for including the corrections from the gauge and 
Higgs sectors. The loop integrals are standard, and the only work is to determine 
the coupling between the particle in the loop and the Higgs bosons J and J7j_ in 
Eq. (^). This approach is easier than including the field dependent masses in the 
formal expression in Eq. (pSf ) and then numerically finding the potential minimum. 
On the other hand, calculating tadpoles alone determines only the first derivatives of 
the one-loop Higgs potential, and does not yield by itself the Higgs potential away 
from the minimum. Fortunately, the minimization conditions are all one needs for 
many analyses. 
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It is crucial to include the one-loop corrections in the effective potential in de- 
termining the vevs. As shown by Gamberini, Ridolfi, and Zwirner^, the tree-level 
Higgs vevs Vi and V2 are very sensitive to the scale at which the renormalization 
group equations are evaluated. Thus it is necessary to determine the proper scale 
at which there are no large logarithms so that the tree-level results are reliable. As 
is well known, there is a simple hierarchy of scales in these theories. As the soft- 
supersymmetry breaking parameters are evolved down from the high scale, the Higgs 
potential evolves so that an asymmetric minimum develops at some scale /iQ. This 
scale is determined by the condition 

ml{^io)ml{^i^) - 5V(Aio) = . (29) 

For Q > /io, the vevs Vi and V2 vanish. For Q < the vevs become nonzero. In the 
supergravity theories under consideration m'jj^ becomes negative, allowing Eq. 
to be satisfied. Figure 3 describes the potential in the regions of interest PSI . 
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Fig. 3. (a) The vevs vanish for Q > hq. (b) For Q « hq, the vevs become nonvan- 
ishing but small, (c) For some scale Q in the range Qo < Q < Ho the vevs have the 
correct magnitudes to give correct electrowcak symmetry breaking, (d) For Q < Qq 
the potential becomes unbounded from below. 



At some lower scale Qq < i^q, the Higgs potential becomes unbounded from below. 
The scale Qq at which this occurs is determined by the condition 

mliQo) + mliQo) - 2\B{Qo)ii{Qo) | = . (30) 

This implies that in the tree-level potential the vevs Vi and V2 must be driven off 
to infinity because the potential becomes unbounded from below. Because the vevs 
evolve from zero at or above the scale /xq all the way to infinity at the scale Qo, the 
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vevs are very sensitive to the scale at which they are evaluated. 

The solution to this conundrum was provided in Ref . . The inclusion of the one- 
loop contributions to the Higgs potential stabilizes the vevs with respect to the scale Q 
at which the parameters (which evolve according to renormalization group equations) 
are evaluated. The standard three cases considered are (i) Msusy < Qo < /^o, (h) 
Qq < Msusy < /^O; and (iii) Qo < /^o < Msusy- In case (i) the scale Qq is determined 



by dimensional transmutation in the sense of Coleman and Weinbergp9[. It was 
initially realized that the one-loop contributions were important in this case, because 
the minimum of the Higgs potential is driven to the flat direction ( "D-flat" ) at tan j3 = 



1 1|30|| , and it was crucial to include the one-loop contribution to lift this degeneracy. 
This yields a light Higgs boson at tree-level (exactly zero mass if tan (3 = 1), which is 
still acceptable experimentally when the one-loop corrections to the Higgs boson mass 



are included|2^. However the predicted SUSY mass spectrum is light and already 
experimentally excluded. Case (ii) has been the subject of much recent work. Case 
(iii) is not of interest since electroweak symmetry breaking does not occur. 

To determine the minimum of the potential, we include the one-loop tadpole 
contributions 

Ti + ATi = , (31a) 
T2 + AT2 = 0. (31b) 

The contributions ATi and AT2 are given in the Appendix. 
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4 Absence of fine-tuning 

The requirement that the supergravity model not be fine-tuned has been recently 
apphed to limit the region of parameter space. This constraint requires that the 
scale of supersymmetry breaking not be too high. Obtaining reasonable criteria for 
declaring a particular theory unnaturally fine-tuned remains a subject of debate. 

The fine-tuning constraint becomes particularly restrictive in the small and large 
tan/3 regions. For small tan/5 (near one), the Higgs potential has its minimum near 
the D-fiat direction. This implies naturally large vevs. Then there must be a can- 
cellation between the two large terms on the right hand side of Eq. (|2^) to obtain 
the experimentally observed value for Mz- Hence for tan/5 — > 1, the supersymmetric 
Higgs mass parameter /i must be tuned ever more precisely - the fine-tuning problem. 
In this section we discuss the various attempts to quantify this constraint. 

The kinds of criteria advocated by other authors are as follows 



Barbieri and Giudice[pT|] introduced a naturalness criteria 

a,- dMl 



< A , (32) 



M| dai 

for various fundamental parameters = mo, mi, /x*"^. A'-' , B'^ to obtain an 
upper bound on the supersymmetric particle masses. They required that A < 
10, i.e. no cancellations greater than an order of magnitude. 

Lopez, et al. define several fine-tuning coefficients e.g. 

5Mz 6fi 
Mz /X 
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They show that a reasonable upper bound on the simplest coefficient implies 
an upper bound on jj,. 

Arnowitt and Nath[0 require that mg < 1 TeV, a condition that is easily applied 
phenomenologically. 



Ross- Roberts]^ and de Carlos and Casas|^ consider the fine-tuning of Mz in 
terms of 

where c is required to be less than some small number e.g. c < 10. Ross 
and Roberts who work strictly with the tree-level Higgs potential argue that 
tan/3 > 2, while de Carlos and Casas argue that the one- loop corrections to the 
Higgs potential ameliorate the fine-tuning. 



Olechowski and Pokorski||10|| look at a full set of derivatives as in Eq. (|33D, ( p4| 



where the Qj are the electroweak scale parameters At, A^, f, tan/5, M^, Mq, 
Mf/, and the Pj are the GUT scale parameters Af, A^, mi, mg, /x'^, , . 
They also find that small tan (3 tends to be more unnatural, and moreover for 
large values of tan /3, near where the top and bottom quark couplings are equal, 
that the model becomes rapidly more fine-tuned as tan/3 is increased. These 
constraints are clearly quite involved. While they test a panoply of fine-tuning 
relations, we feel they are overly complex for such a qualitative and arbitrary 
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notion as naturalness. Therefore we abandon this notion in favor of a more 
intuitive definition similar to Lopez, et al.[3. 



• Castano, Piard, and Ramond|Tl| choose a numerical definition in which the 
number of iterations the computer has to find a solution is limited. It is not 
obvious how this algorithm compares quantitatively to those defined above. 

Our physical definition of naturalness is simply \fi{Mz)\ — \fi{mt)\ < 500 GeV. 
A measure of the reasonableness of this definition is the effect that small changes 
in /i have on Mz- From the tree level equation for (see Eqn. (|2^) ) it is readily 
apparent that larger values of become more unnatural. In Figure 4 we plot the 
dependence of on fi and B for both the tree-level calculation and the full one-loop 
contributions to the Higgs potential. It can be seen that the one-loop contributions 
reduce the fine-tuning to an extent. This comparison can be related to the plots of 
the vevs as a function of scale Q as discussed in Ref. P2 . 
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-230 -220 -210 -200 -190 -180 -170 460 -450 -440 -430 -420 -410 -400 

B(Mgut) B(Mgut) 
Fig. 4. The change in Mz with /i and B for the fixed-point solution mt{mt) = 160, 
tan (3 = 1.47, niQ = 0, = 0, and two different values of mi. The values of ^{Mq) 
meet our naturalness criterion |/x(M^)| < 500 GeV. The solid (dashed) curves are 
the results at one-loop (tree) level. The case < gives comparable curves. 



The plots in Figure 4 correspond to a low-tan/? fixed-point solution[15|-|18 



m 



such cases the tree-level and one-loop-level values turn out to be comparable for 
either /i and B in the region defined by our naturalness criterion (only the degree of 
fine-tuning changes). Consequently, including the one- loop corrections in the Higgs 
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potential does not have a critical impact on the phenomenology. This result does not 
extend to other regions of the m^-tan j3 plane since there our criterion for naturalness 
implies a larger allowed range of mo and mi in which the one-loop contributions can 
change and B significantly (see section 7). 

One can consider quantitative fine-tuning criteria analogous to those considered 
above, 



5Mz 




Mz 


~ '''' I^G 


5Mz 




Mz 





(36) 



with e.g. |cj^|, \cb\ ^ 30, where the derivatives on the left hand side are obtained at 
the physical Z mass scale and the derivatives on the right hand side are at the GUT 
scale (denoted by the G superscript). Since the RGE equation for n is proportional 
to /I, the value of S/i/ii is scale independent, but SB/B depends on scale. Table 1 
gives the values of |c^| and \cb\ determined for the tree- level and one- loop curves for 
the low-tan P fixed point solution of Figure 4. 



mi 


|cmI, (/^ > 0) 


|cb|, {h > 0) 


|c^|, in < 0) 


\cb\, (/U < 0) 


100 GeV (loop) 


8.8 


8.2 


7.3 


5.2 


100 GeV (tree) 


13.5 


10.9 


13.5 


6.3 


200 GeV (loop) 


25.6 


27.0 


20.3 


16.9 


200 GeV (tree) 


57.0 


46.8 


57.0 


27.6 



Table 1: Values of |c^| and \cb\ obtained at tree 
and one- loop levels for the fixed-point solution of Figure 5. 
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Note that inclusion of the full one-loop contribution substantially reduces the fine- 
tuning constants |c^| and \cb\- Our entries for |c^| are somewhat larger than those 
found in Ref. because our model has a value of tan/3 that is closer to tan/5 = 1. 

5 Models 

The introduction of supersymmetry introduces many new unknown parameters to 
the standard model. The advantage of the popular supergravity models is that this 
number of new parameters is reduced to five or less. The models discussed here 
should only be viewed as examples of possible supersymmetry breaking scenarios. 
Some features may be more general, however. 

A. General model[j5[] 

The universal parameters at the GUT scale are tUq, A'^, mi, /i*^, B'^. In the 
minimal supergravity model, these five parameters describe the higgsino and gaugino 
sectors. The universality of the scalar masses at the GUT scale provides for the 
suppression of dangerous fiavor changing neutral currents involving the squarks of 
the first two generations. 

B. No-scale 

In no-scale models two of the five parameters are zero at the unification scale, 

mo = 0, = . (37) 

Thus the scalar fields are massless there, and mi is the sole origin of supersymmetry 
breaking. 



C. Strict no-scale 1 30, 33 
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The strict no-scale model is a version of the no-scale model with 

= 0, (38) 

at the unification scale. 



D. Dilaton 20 



When the dilaton S receives a vev, one encounters a breaking of supersymmetry 
that is of a different nature than that of the minimal supergravity scenarios described 
above. The dilaton F-term scenario leads to simple boundary conditions for the 
soft-supersymmetry parameters 

TJin = —=mi, = —mi . (39) 

This model therefore has only three parameters. When it is required that fi receive 
contributions from supergravity only, the additional unification constraint 

B^ = 2mo , (40) 

is obtained. The dilaton version of supersymmetry breaking has been studied in the 
MSSM in Ref. and for the fiipped SU(5) model in Ref. [g§ . 



E. String-Inspired 

Supersymmetry breaking in strings is a nonperturbative effect, since supersym- 
metry is preserved order by order in perturbation theory. Very little is known about 



nonperturbative effects in string theory. Recently the authors of Ref. have pro- 
posed to parametrize our ignorance of the exact nature of the breakdown of super- 
symmetry. The dilaton breaking scenario above is a specific case of more general 
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scenario of supersymmetry breaking in which the moduh fields also receive a vev. 
If one restricts oneself to the case where only one T field and the dilaton S get vevs, 
then the amount of SUSY breaking that arises from each sector can be parametrized 
by the "goldstino angle|2T|" 6. The dilaton breaking case corresponds to sin^ = 1. 



The angle 6 is constrained by low-energy phenomenology since purely dilaton break- 
ing gives a universal boundary condition for the scalar masses, and the breaking of 
supersymmetry when the moduli field gets a vev will gives rise to FCNCs in the 



low-energy theory. According to Ref. the more general case, where substantial 



contributions to supersymmetry breaking arise from the moduli field getting a vev, 
is not ruled out. 

The unification scale in the string-inspired model is roughly an order of magnitude 
higher than the scale at which the gauge couplings unify in the MSSM. Presumably 
large threshold corrections due to non-degenerate GUT particles could account for 
this discrepancy. 



F. Large tan/3 scenario [|36[] 



The correct electroweak symmetry breaking does not occur for too large values 
of tan/3. If tan/3 ^ mt{mt)/mf,{mt), then the bottom quark Yukawa drives the 
Higgs masses parameter ml negative first (instead of m| from the top quark Yukawa 
coupling). For tan/3 close to this limit considerable fine-tuning is required to get the 
correct electroweak scale. This situation is ameliorated somewhat with the inclusion 
of the one-loop corrections in the effective potential [p!0[] . 
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6 Ambidextrous Approach to RGE Integration 

Previous RGE studies of the supersymmetric particle spectrum have evolved from in- 
puts at the GUT scale (the top-down method |^ |rT|) or from inputs at the electroweak 
scale (the bottom-up approach |TD|). Our approach incorporates some boundary con- 
ditions at both electroweak and GUT scales, which we call the ambidextrous approach. 
We specify rrit and tan/3 at the electroweak scale (along with Mz and Mw) and mi, 
rriQ, and A'^ at the GUT scale. The soft supersymmetry breaking parameters are 
evolved from the GUT scale to the electroweak scale and then fi{Mz) and B{Mz) 
(or nirrit) and B{mt)) are determined by the tadpole equations at one-loop order. 
Subsequently /i and B can be RGE-evolved up to the GUT scale. This strategy is 
effective because the RGBs for the soft-supersymmetry breaking parameters (see the 
appendix) do not depend on /i and B. This method has two powerful advantages: 
First, any point in the rrit - tan/3 plane can be readily investigated in specific super- 
gravity models since rrit and tan j3 are taken as inputs. Second, the tadpole equations 
Eq. (|79a| - [79b|) are easy to solve in the ambidextrous approach. The Ti equation can 
be solved iteratively for fi{Mz), and then the T2 equation explicitly gives B{Mz)- We 
stress the numerical simplicity: no derivatives need be calculated and no functions 
need to be numerically minimized. 

We now describe our numerical approach in more detail. Starting with our low- 
energy choices for rrit, tan/?, 0:3, and rrif, (and using the experimentally determined 
values for cti, 02 and mT-[^]), we integrate the MSSM RGEs from rrit to Mq with 
Mq taken to be the scale Q at which ai{Q) = a2{Q)- We then specify mi, mo, and 
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A*^, and integrate back down to where we solve the tadpole equations for nirrit) 
and B{mt). We can then integrate the RGEs back to Mq to obtain fi^ and B'^ at 
Mq. a few remarks are pertinent: 

1) We integrate the two-loop MSSM RGEs for the gauge and Yukawa couplings |T3|] , 
but only the one-loop MSSM RGEs (as given in the appendix) for the other super- 
symmetric parameters. We retain the important two-loop gauge and Yukawa effects 



(only the two-loop gaugino RGEs exist |^ and we desire to be consistent with regard 
to the order for the soft-supersymmetry breaking parameters). 

2) Since we are working only to one-loop order in most RGEs, we neglect threshold 



effects at both the GUT||3^-|^ and supersymmetry scales. To properly take into 
account threshold effects requires that the appropriate beta functions be calculated 
at every (s)particle threshold, and since the complete supersymmetric mass spectrum 
and the GUT particle spectra are generally not known a priori, the calculation and 
use of these beta functions is extremely daunting. In any event, it should not be 
critical to incorporate these threshold effects into a one-loop calculation since two- 
loop contributions will be comparable to threshold effects. 

3) We take the lower bound of our integration at rrit instead of Mz for several 
reasons. As shown by several groups[^, ^, |10], |lT|, inclusion of the one-loop effects 
into the effective potential makes electroweak symmetry breaking roughly indepen- 
dent of scale; the scale Q = rrit is roughly the value for which the large logs cancel 
among themselves in the one-loop corrections to the minimization conditions. We 
choose rrit as the boundary since the RGEs (in particular for the gauge and Yukawa 
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couplings) are simple at scales above m^, and it is non-trivial to extend them below 
rrit. In addition, the choice of rrit facilitates comparison with previous work on gauge 
and Yukawa unification and fixed points. 

7 Results 

We discuss the supersymmetric spectrum and phenomenology for several represen- 
tative points in the rrit ^ tan/5 plane. For the most part we focus on the low-tan /3 
fixed-point region since it is very attractive to explain the large top quark mass as 



a fixed point phenomenon |T^-[|T^. Moreover, the supersjTiimetric spectrum in this 
region is largely unexplored, probably due to fears of excessive fine-tuning. How- 
ever, as addressed in Section 4, these fears are not necessarily justified; there remains 
substantial viable parameter space for which fine-tuning does not pose great con- 
cern, particularly with the inclusion of the full one-loop corrections to the effective 
potential |I0|] . 

The rut - tan/5 parameter space can be divided into several distinct regions, as 
shown in Figure 5. 



28 



60 



50 



40 



^ 30 



20 



10 







100 



M 



susY = "^ ; a3(Mz) = 0.12 



(e) 



m (m) =4.25 GeV 
b b 




120 



140 



160 



180 



200 



m^(m^ GeV 

Fig. 5. The allowed mt - tan f3 parameter space assuming Yukawa unification 
Xb{MG) = At-(Mg'). [^]. The shaded area indicates the region for which mb{mh) = 
4.25ib0.15 GeV. Points representative of distinct regions within this parameter space 
are denoted with labels (a)-(e). 

We discuss the supersymmetric mass spectrum for each of these regions. Unless 
otherwise specified, we take A"^ = 0, a3{Mz) = 0.120, and mi,{mh) = 4.25 GeV. The 
qualitative behaviour in each region should not depend greatly on these parameters. 



(a) Low-tan /3 Fixed Point 

As a typical example of the low tan f3 fixed values region we consider the point 
mt{mt) = 160 GeV, tan/3 = 1.47 (for which XtiMc) = 2.7). We aim to deter- 
mine the GUT-scale parameter space for which this solution can be obtained 
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from the minimization of the effective potentiaL Using the tadpole method, we 
explore a grid of mg and mi values and apply both experimental and natural- 
ness bounds. For the lower experimental limits, we adopt the values listed in 
Table 2 following Ref. g|. 



Particle 


Experimental Limit (GeV) 


gluino 


120 


squark, slepton 


45 


chargino 


45 


neutralino 


20 


light higgs 


60 



Table 2: Approximate experimental bounds that we apply in Figure 6. 



Figure 6 shows the allowed parameter space for both signs of /x along with 
the most restrictive constraints in each case. The contours of constant |yu| are 
ellipses in the mo — mi plane for \n\ » Mz- 
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Fig. 6. The allowed mo, m-i region is shaded for the low-tan /3 fixed point mAmt) = 

2 

160 GeV, tan/3 = 1.47 solution with (a) /i > and (b) fi < 0. The experimental 
bounds in Table 2 and the naturalness bound \fi{mt)\ < 500 GeV are imposed with 
= GeV. A semi-quantitative dark matter constraint (given by Eq. (]45|)) is also 
shown. 
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Note that the /i < case has more available parameter space; it is also slightly 
more natural, as indicated from the fine-tuning constants given in Table 1. 

We have indicated the light scalar higgs experimental limit with a dashed 
line in Figure 6; care must be taken when enforcing this particular constraint 
since the allowed parameter space is somewhat sensitive to the exact limit. 
Moreover, the rrif^ bound includes only the one-loop quark-squark contributions 



given in reference and it is expected that inclusion of the chargino and 
neutralino contributions can affect the mass of the light scalar higgs by a few 
GeV[^]. In general an increase in rrih (as in all the experimental constraints) 
tends to increasingly restrict the allowed parameter space, while a decrease in 
m/j has the opposite effect. 

Overall we find substantial phenomenologically viable parameter space, 
especially for /i < 0. However the maximal values of the GUT parameters mo 
and mi are not large (mo ^ 350 GeV, mi ^ 225 GeV) implying a rather light 
low energy supersymmetric mass spectrum. Also included in Figure 6 is the 
semi-quantitative dark matter constraint of Drees and Nojiri (see Eq. (^5]) 
below). For this low tan/3 fixed-point case it implies that mo ^ 250 GeV, 
though this approximate bound ought not to be taken strictly. 

We now investigate the supersymmetric particle mass spectra dependence 
on mo and mi independently for this low tan/3 fixed-point solution. Fig- 
ure 7 shows the dependence of the supersymmetric spectrum on mi in the 
noscale model (mo = 0). For the squarks and sleptons we plot the light- 
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est mass eigenstates; in addition we plot the heaviest stop, m^^, for refer- 
ence. We label the chargino and neutralino masses such that M± < M± 
and M^o < M^o < M^o < M^o. 
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Fig. 7. The low tan /3 fixed point solutions for (a) /i > and (b) /i < with 
mo = GeV and AP = 0. The experimentally excluded region includes all experi- 
mental constraints except for the bound on m/i, since it is sensitive to chargino and 
neutralino contributions [H^l. 
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Figure 8 shows all the squark and slepton masses for the same low-tan (3 
fixed-point solution with > 0. Note that the squarks of the first two gener- 
ations can be heavier than those of the third; the up and charm squarks are 
degenerate as are the down and strange squarks. The slepton masses are ap- 
proximately generation independent in this case, though this need not be true 
in general (e.g. see Table 3 below). 
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Fig. 8. The squark and slepton masses for the low-tan [3 fixed-point solution in the 
no-scale model with /i > 0. The solid (dashed) lines correspond to the third (first 
and second) generation. The excluded regions are the same as in the previous figure. 



Figure 9 illustrates the dependence of the supersymmetric spectrum on 
niQ. (Here we take mi = 150 GeV.) The mass of most of the SUSY particles 
increase with increasing mg (see, e.g. Eq. (P)). 
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Fig. 9. The low tan/3 fixed point solutions for (a) > and (b) < with 
mi = 150 GeV and A*^ = 0. The shaded area denotes the region excluded by our 
naturalness criterion. 
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We also give qualitative descriptions of the allowed parameter space in the 
other significant rrit - tan /? regions. 

• (b) Medium-tan /3 Fixed Point 

The allowed mo ^ '"^i parameter space is substantially larger in this case than it 
is for the low-tan /3 fixed point. Our naturalness condition allows substantially 
larger values of mo and mi (mo ^ 725 GeV and mi < 325 GeV for mt{mt) ~ 
192 GeV, and tan/3 = 15); however, dark matter constraints will still require 
mo ^ 300 GeV. For /i > 0, experimental bounds on nif^, and m^^^^. push 
the lower bound for mo and mi up slightly. For ^ < 0, experimental bounds 
for mo and m^^ also become more restrictive, but the m^^j_ and m^^ constraints 
become less restrictive. In both cases the constraint from the lightest scalar 
higgs mass becomes less restrictive; even in the /i > case, it will not play 
an important role in limiting the allowed mo and mi region. To summarize, 
the medium-tan /3 fixed point region allows larger values of mo and mi from 
our naturalness constraint while the experimental restrictions on this parameter 
space do not change much from the low-tan P case. 

• (c) High-tan /3 Fixed Point 

This region describes the SO (10) relation = A;, = A^ where Aj ^ 1. There is 
not much parameter space remaining without weakening our naturalness con- 
dition. For the case mt{mt) — 178 GeV, tan/3 = 61 (with mo = 400 GeV, 
mi = 400 GeV, and ^x{rnt) ~ 575 GeV) the particle spectrum is given in the 
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following table. 



Particle 


Mass (GeV) 


gluino 


1078 


stop, sbottom 


751,900; 763,881 


up squarks, down squarks 


1029,1060; 1026,1063 


stau, tau sneutrino 


183,454; 417 


other sleptons 


431,494; 487 


charginos 


323,590 


neutralinos 


167,323,579,588 


higgs: mA,mH±,mH,mh 


364,377,363,131 



Table 3: Particle spectrum for rrit = 178 GeV, tan/3 = 61 
(where mo = 400 GeV, mi = 400 GeV, ^ 0). 



As before the above particle spectrum is calculated at the scale m^. We obtain 
no natural solutions for < 0. 

• (d) Low-tan /3, not a Fixed Point 

This region has a large amount of viable parameter space; naturalness bounds 
allow substantially higher values of mi (mi ^ 330 GeV for mt{mt) = 160 
GeV with tan (3 — 3), and experimental constraints do not further restrict this 
parameter space to any great extent (though the sneutrino and chargino bounds 
are pushed upward somewhat). In fact, the higgs constraint is weakened a great 
deal for > 0, allowing relatively low values for mo and mi. Moreover, the 
light stop constraint is less important in the /i < case. 
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• (e) High-tan/?, not a Fixed Point 

The allowed parameter space is reduced by the lightest stau constraint (which 
cannot be the LSP), though some parameter region remains. The allowed pa- 
rameter space is bounded by chargino, stau, dark matter, and our naturalness 
constraint which give 180 ^ mo ^ 300 GeV and 85 ^ mi ^ 400 GeV for 

2 

mj(mf) = 160 GeV, tan/3 = 45. The light higgs and the light stop constraints 
are not important for either sign of fi. 

In addition, we varied from —500 to +500 GeV in the low tan P fixed point 
case; we found little change in the resulting parameter space except that the light 
stop constraint is more (less) restrictive for A'-' negative (positive) and yU < 0. The 
fixed point solution in radiative electroweak symmetry breaking has also been studied 



recently in Ref. [113 



A critical constraint [^, Q on the supersymmetric spectrum is the rare decay 
b ^ s'-f. We remark here that regions of the parameter space illustrated in the 
previous figures are not ruled out by this constraint. This will be the subject of a 
forthcoming paper [^. 



8 SUSY Mass Spectrum Correlations 

For smaller values of tan/3 it is clear from the tree-level expression Eq. (^) that 
is usually large compared to the the electroweak scale Mz- Furthermore the fine- 
tuning problem in this situation is softened when the one-loop contributions to the 
Higgs potential are included. For values of /i just a few times larger than Mz, the 
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particle spectrum is governed by certain asymptotic behaviors which we discuss in 
this section. 

As discussed previously, the gaugino masses are related (through one-loop order) 
by the same ratios that describe the gauge couplings at the electroweak scale. This 
observation, together with the fact that is large, yields simple correlations between 
the lightest chargino and neutralinos and the gluino0, |5^, namely 





~ Ml, 






(41a) 


^x? 


~ M^o 

X2 


~ M2 = 


—Ml ~ 2Mi ~ M^o , 

ai ^'^ 


(41b) 


m-g 


= Ms 


= ^M, = 

^2 




(41c) 



where the quantities in these equations are evaluated at scale nit. The heaviest 
chargino and the two heaviest neutralino states are primarily Higgsino with 

M^± ~ M^o ~ M^o ~ . (42) 

(43) 

The lightest Higgs h has small mass for tan j3 near one at tree-level by virtue of the 
D-flat direction; its mass comes from radiative corrections p3, The heavy Higgs 



states are (approximately) degenerate ~ Ma because at tree- level Ma = ~ 
—Bfi is large. 

The squark and slepton masses also display simple asymptotic behavior at large 
The first and second squark generations are approximately degenerate (though 
not degenerate enough to ignore their contributions to the minimization of the effec- 
tive potential). The squark and slepton mass spectra are shown in Figures 7 through 
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9. The splitting of the stop quark masses grows as \fi\ increases. This sphtting of the 
sbottom states does not change much with fi for small tan f3. 

9 Dark Matter 

The neutralino as the LSP is an ideal candidate for the dark matter since it is sta- 
ble and interacts weakly. The MSSM utilizes R-parity conservation so the lightest 
neutralino must annihilate to ordinary matter {xx — ^ R — even matter) to a suffi- 
cient extent to avoid overdosing the universe[^]. For a bino-like LSP, dark matter 



considerations put an upper bound on the parameter ttiq. We adopt the conservative 
viewpoint that the contribution of the LSP alone to the dark matter of the universe 
must be less than the closure density. Roberts and Roszkowski [^] apply an addi- 
tional constraint in which the neutralino is required to make up a substantial fraction 
of the dark matter; this requirement provides a lower bound on tuq as well. The 
recent results from COBE suggest that the dark matter is a mixture of hot and cold 
dark matter. Although it may be simpler to assume that all of the cold dark matter 
is composed of one contribution, it is perhaps premature to assume this. We remark 
that the recent exciting claims of experimental evidence for dark matter in nearby 
galaxies |]5^, ^ only solves the local baryonic dark matter problem[Q. The origin of 



the nonbaryonic dark matter needed to close the universe is still unknown. 

The typical situation in the low tan/3 fixed point solutions is that \fi\ » and 
consequently the lightest neutralino (which is the LSP) is predominantly gaugino; 
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indeed the LSP is predominantly bino. The neutrahno mass matrix is 



/ Ml 


—Mz cos P sin 9, 
\ Mz sin P sin 





M2 
cos P cos 
— M^sin/^cos ^. 



—Mz cos /? sin 6* 
Mz cos /3 cos 



w 



w 



Mz sin /3 sin 6^r 
—Mz sin /3 cos 9^, 





w 







(44) 



For >> M2 the hghtest two neutrahnos are predominantly bino and wino, and 
hence the bino and gaugino purities are high. In this case any mass limit on the 
lightest neutrahno from Z decays at LEP disappears since the Z couples only to the 
Higgsino component of the neutrahno. Figure 10 gives the bino and gaugino purities 
for the low-tan P fixed point solution in the no-scale model, corresponding to Fig. 7. 
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Fig. 10. The bino and gaugino purities for the low tan/3 fixed point solution with 
mt = 160 GeV, tan/3 = 1.47 in the no-scale model with (a) > and (b) n < 0. 
Shaded regions are forbidden by experimental and fine-tuning considerations. 
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Given that the solutions are comfortably in the high bino purity region, we apply 



to obtain the line corresponding to Qh'^ = 1 in Figure 6. This formula is based 
on the observation that for the bino-like LSP the annihilation rate is dominated by 
the sleptons, and it neglects a possible enhanced annihilation rate that may occur if 
there are significant s-channel pole contributions. The bino and gaugino purities for 
non-zero mo (in particular the dilaton model) are similar to the above figures. 

10 Conclusion 

The motivation of this work has been to distill the interesting supersymmetric phe- 
nomenology of the low-tan (3 fixed-point region that can explain the origin of a large 
top quark mass. The RGEs are solved with some boundary conditions taken from 
both GUT and low energy scales. The minimization conditions on the effective po- 
tential are obtained with the tadpole method. 

Our principle findings can be summarized as follows: 

• Solutions with a At fixed point, rrit ^ 170 GeV and radiative breaking of the 
electroweak symmetry breaking are allowed. These solutions are characterized 
by relatively large values of the supersymmetric Higgs mass parameter which 
implies that the supersymmetric particle spectrum displays a simple asymptotic 



the semi-quantitative constraint of Drees and Nojiri [^ (valid roughly for > mi 
M^o > 60 GeV) 
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behavior. The solutions also meet the naturalness criterion \fj,{Mz)\ < 500 GeV 
for both signs of /i. 

Representative sparticle mass spectra are presented for the At fixed point solu- 
tions. 

Over most of the GUT parameter space for the low-tan /3 fixed-point, the gaug- 
ino masses exhibit show simple correlations due to the relatively large value 
compared to M2. The heaviest chargino and the two heaviest neutrahnos 
have masses approximately the lightest chargino and the second lightest 
neutralino have masses approximately M2; the lightest neutralino (LSP) has a 
mass approximately Mi ~ M2/2. The lightest Higgs obtains its mass almost 
entirely from radiative corrections, and the states H, H^, A are relatively heavy 
and approximately degenerate. 

In the early universe the LSP will annihilate sufficiently neglecting s-channel 
pole annihilation for most of the parameter space (mo ^ 300 GeV) so as not to 
overdose the universe. 

The values of // and B derived from the one-loop Higgs potential analyses are 
very similar to the tree- level results in the low-tan /? fixed-point region when the 
parameters Mz and tan (3 are taken as input. However, the one-loop corrections 
to the Higgs potential somewhat ameliorate the fine-tuning problem. 

The tadpole method is a convenient way to calculate the one-loop minimization 
conditions. We have obtained these conditions in an analytic form including all 
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contributions from the gauge-Higgs sector and matter multiplets. 
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12 Appendix 



Renormalization Group Equations 

The renormalization equations for the gauge couplings and the Yukawa couplings 
to two-loop order can be found in Ref. |T^. In the most general case the evolution 
equations involve matrices. For example the Yukawa couplings form three-by-three 
Yukawa matrices: U for the up-type quarks, D for the down-type quarks, and E for 
the charged leptons. Similarly the soft-supersymmetry breaking parameters form the 
matrices Mq^, M^/^ , Mj^^, and M^;^. Finally there are in general matrices for 

the trilinear soft-supersymmetry breaking "A-terms": Au, Ad, and Ae- It turns out 
to be useful to define the combinations Uajj = AuijUjj, etc. in the matrix version 
of the RGB's. Then the evolution of the soft-supersymmetry parameters (with our 
convention for signs) is given by the following renormalization group equations [H^ 



~dt 



dEj, 

dt 



dMi 


2 


dt 


167r2 


dVA 


1 


dt 


167r2 




HP 



(46a) 



-2(DADtU) + (DDtUA)} , (46b) 

+ {[4(DaD"^D) + 6Tr(DAD"f)D] + [^{T>T>^T>^) + 3Tr(DD"^)DA" 
+2(UAUtD) + (UU^Da) + 2Tr(EAEt)D + Tr(EEt)DA} , (46c) 



167r 



167r2 



?>gi + 3(72' Ea + 2 3(7i'Mi + ?>giM2]E 
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+{ [4(EaE^E) + 2Tr(EAE^)E] + [5(EE^Ea) + Tr(EE^)EA 
+6(DADtE) + 3(DDtEA)} , (46d) 
^ : (IglM, + 3glM2 + Tr(3UUA + 3DDa + EEa)) , (46e) 



167r2V5' 

^ ^ - 3(72' + Tr(3UU^ + 3DDt + EEt)) , (46f) 



167r2 V 5 



+3Tr(D(M|^ + M^jDt + M|^DDt + DaD^) 

+Tr(E(Mi^ + MljEt + M^^EEt + EaEJ,)) , (46g) 



167r2 V 5^ 

+3Tr(U(M|^ + M^^jUt + M|^UUt + UaU^)) , (46h) 

T^,{ - ^gfM'i - SgM - ^-^alMi 



+i[UUtM^^ + Mj^UUt + 2{VMl^V^ + m^^UUt + U^U^)] 
+^[DDtMj^ + M|^DDt + 2(DM^^Dt + m^DDt + D^dI)]) , (46i) 

+ [V^lJMl^ + M^^UtU + 2(UtMj^U + m^^UtU + U^U^)]) , (46j) 
167r2V 15^^ ^ 3 

+[DtDM^^ + M^^DtD + 2(DtM|^D + mjj^B^B + 0^0^)]) ,(46k) 



r^,{-lgM-3gM 

+ ^[EEtMi^ + Mi^EEt + 2(EM|^Et + m^^^EEt + E^E^)]) , (461) 
+ [E^EMl^ + Ml^E^E + 2(EtMi^E + m^^E^E + E^E^)]) , (46m) 
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For our purposes it is sufficient to consider these equations keeping only the leading 
terms in the mass hierarchy in the three generation MSSM. The resulting renormal- 



ization group equations]^ are given below to leading order. 



dt 
dAt 

dt 
dAb 
IT 
dAr 
IT 

dB 

m 

dfi 

It 
dMj^ 

dt 
dt 
dt 

dt 
dMl 

dt 
dMl 

dt 
dM^ 
dt 



2 

167r2 


bi9-Mi , 


(47a) 


2 


{Y.Cig^M, + 6\',A + XlA,) , 


(47b) 


2 

167r2 


( ^ 4g^M, + GXlAt, + A + A.) , 


(47c) 


2 

167r2 


{J24g^M, + 3XlA, + AXlA,) , 


(47d) 


2 


(^glM, + 3glM2 + 3XlA, + 3A?A, + A^A.) , 


(47e) 


fi 
167r2 


{-lal- + 3A? + 3X1 + A^) , 


(47f) 


2 

167r2 


( - ^iMl - 3glMl + 3XlX, + A^X.) , 


(47g) 


2 

167r2 


( - IglMl - 3^72^M| + 3A?X,) , 


(47h) 



167r2 
2 

167r2 
2 

167r2 
2 

167r2 
2 

167r2 



1 16 

-glMl - 3glMl - -glM^ + XlX^ + A^X, 
f/.Ml - ^^glMl + 2XlX,) , 
y/^MI - ^-^glMl + 2A^X, 



glMl - 3glMi + A^X. 



2 T\l2 



^^glMl + 2A^X. 
5 



and for the two light generations 
dAu 



dt 
dAd 

dt 
dAe 

~dr 



2 
2 
2 



(47i) 
(47j) 

(47k) 
(471) 

(47m) 



{j:^':9iM.+xiA+\xiA^) 



(48a) 
(48b) 
(48c) 



50 



where 



■ ( - \m - ZglMl - ^glMl) , (48d) 



dt 167r2V 15' 

- {- T^glMl - -glMi) , (48f) 



167r2 V 15 

dM? 2 / 3 



dt 167r2 V 5 



( - -glMl - ZglMl) , (48g) 



^ML ^ 2 . 12 
167r2 V 5 



^ - — (48h) 





= (f,l,-3), 




(49a) 


Ci 


" 45'^' 3^' 




(49b) 




" 45'^' 3^' 




(49c) 


4 


= (^3,0), 




(49d) 


Xt 






(49e) 








(49f) 


Xr 






(49g) 



Here the factors q, c^, and c'l are given by a sum over the fields in the relevant Yukawa 
coupling, e.g. q = J2fCi{f) = ^(-^^2) + CiiQ) + CiiU^)- The coefficients in front of the 
gauge coupling parts of Eq. (41)-(43) can be understood from the quantum numbers. 
For the fundamental representations of SU{N) there is a factor of {N"^ — 1)/N and 
for the hypercharge U{1) one has -^Y"^ (with hypercharge suitably normalized, e.g. 
n« = 2). 

One-Loop Effective Potential 

51 



We summarize here the tools needed to construct the one-loop minimization con- 
ditions. The necessary ingredients are the field dependent particle masses; since we 
are calculating the tadpole diagrams, we need the particle masses at the potential 



minimum and the Higgs couplings. The tadpoles are calculated in the DR renormal- 
ization scheme [p5|] . 

We present here the contribution from the third generation (s) particles; the contri- 
butions from the other generations can be obtained with obvious substitutions. The 
top and bottom squark and the tau slepton mass matrices (at the potential minimum) 



are 



Mql + + - Ml) cos 2/3 mt{At + /i cot /3) 

mt {At + ^l cot M^^ + - i {M^ - M|) cos 2(3 



(50a) 



M^^ + ml -li2M^ + Ml) COS2P mb(^ + /i tan/3) \ 

mbiAb + fitanp) + + |(M2, - M|) cos2/3 ) ' ^''"^ 



Ml^ +ml- \{2Mlr - M|) cos 2/3 m^(A^ + /itan/3) 

mr {Ar + ^l tan /3) M^^ + ml + {Mir - M|) cos 2/3 



(50c) 



which are diagonalized by orthogonal matrices with mixing angles 9i, 6^, and 6f. The 
mass eigenstate for the massive third generation sneutrino is 



ml = M^^ + cos 2/3 . (51a) 
The relevant Higgs couplings to the squark eigenstates are 

wMz \f cosHt \fl • 2. V / sinHt \, . 2. ; 



cos 9^ 



(52a) 



52 



V{Jf2f2) 



+ 



cos 9 



w 



cos^ 9b 
sin^ 9b 



Cb sin 



tan (3 T sin 26^6 [^6 tan (3 + n] 



2Mw 



igMz 
cos ^ 



2M, 



VI/ 
COS^ 

sin^ 6t 



- + sin ) + 



sin^ 9b 



^^"^■"^ tan /3 T sin 29r [Ar tan/3 + /x] 



sm 
cos^ 



(cft sin^ 61^^) 



(52c) 
(52d) 



{ct sin^ ^ly) 



(52e) 
(52f) 



where e*, Cb, and Cr are the electromagnetic charges 2/3, —1/3, and —1 respectively. Notice that 
the D-terms do not contribute to the coupling of Jj_ to the squarks. The mixed couplings (e.g. 
V{J'iii2)) obviously do not contribute to the tadpole. Calculating the tadpole, and making use of 
the relations 



sin 2^4 
cos 29 1 
sin 26b 
cos 26b 
sin 29r 
cos 26r 



_ 2mt(At+Mcot/3) 



2 

ti 

2 



2 

2 

- TO- 



2 

TO^ 

tl 



2 

TO^ 

12 



2to? 



Ol 



TO? 



b2 
^2 



TO? — TOj 

Ol Or 



TO? 

01 



TO? 

02 



m% 



2 



(M2^ - M^) + i cos 2/3(8M2, - 5M|) 



2 

tl 



t2 



2mb{Ab + /xtan/3) 



TO? 



TO? 



bi 62 

M2 )_icos2/3(4M2,-M|) 



2 2 
m~ — ra% 
61 ^2 



2tot-(j4^ + /utan/3) 



TO? 

Tl 



to: 



{Ml^ - M^J - i cos 2/3(4M2, - 3M|) 



™li 



(53a) 
(53b) 
(53c) 
(53d) 

(53e) 
(53f) 



one arrives at the top and bottom quark-squark contribution to the minimization conditions: 



(9) 




-Ml 



- hMl) \\(Ml^ - Ml) + 1 cos2/3(8M2, - 5M|) 



+2ml{{iicot(3f -4) 



-m 



bi, 



In 



2 

01,2 



- 2ml + 



53 



± 



m? — m? 

bi 62 



i(4M^-M|)fl(M2^ 



M? 

Or 



12 



cos2/3(4M^ - M 



+ 



+2mg((Mtan/?)2- Ag) 
3 



(54a) 



3 f 



mf I In 



tl,2 



In 



2 

m£ 

2 



1 ) cot /3 + ( In 



n2 



1 tan /3 



+ mi 



In 



2 

mt 



g2 

- 1 I tan /3 



1 cot/3 



m? ± 



mi 



-ml{At + /icot/3)(At + /xtan/3) 



— 

Oi b: 



-ml {Ab + cot /3) + /Lt tan /3) 



..(54b) 



Also the tau lepton-slepton and sneutrino contributions are 



(0 



1 



—m^ In 



mz 



+ 



mlMz' 



± 



In 



2 

m% 



1 



2mi 



^^1 
2 ^ 



™ll - ™f2 



(4M^ - 3M|) 



^(M|^-M4)-i cos 2/3(4M^ 



3M|) 



+2m^(^(/itan/3)2 - 



(55a) 



1 



ml± 



m^(ln-^-l)tan/? 



T^l,2 



In- 



1 tan/3 



2 

TO- 



- TO^ (^r + A* cot /3) ( At- + /U tan /3) 



(55b) 



There are similar contributions from the first and second generations. Most of these terms in 
ATi and all of the terms in AT2 are proportional to some powers of the quark or lepton mass, 
which is negligible in the light generations. However, there exist contributions to ATi which are 
proportional to Mf and M^; these are zero only in the limit in which the squarks (or sleptons) 
are degenerate within each generation. This is not necessarily a good approximation; we find that 
the light squark/slepton contribution can be larger than the gauge boson contribution (sec below), 
especially for moderate or large values of TOq and mi. It is therefore important to include the light 
squarks and sleptons in a full one-loop analysis. Explicitly, the light squark and slepton contribution 
is 



(2)3 



— TO 



di ■ 



In 



2 

TO - 
<ii,2 



-A 



iM|±(i)i(8M^-5M|) 



^^I±(^)^(4M2,-M|) 



(56a) 
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AT, 



ill) 



(2)1 



,n?,(l„f -1 



-tij, „ In 



mi 

ei,2 



iM|±(i)(4Af2,-3M|) 



AT,^'"'^ = , AT,^"^ = , 



(56b) 
(56c) 



where the factor of two includes both hght generations. 

If we neglect the contribution from the bottom quark and from the D-term contributions to the 
squark masses and couplings, the equations above reduce to 

f (to? ) — f (to? ) / s 
2/(TO^) - fiml) - /(TO?J + ii-iii_^^((^cot/3)^ - 



AT, 



(9) 



AT. 



3to2 



3toj cot /? 
8n^v 



TO? — TO? 

tl 12 



2f{fnf}-f{ml)-f{ml) 



fiml)-f{ml) 



9 9 
TO- — TO- 
tl 12 



{At + n cot /3) {At + /i tan ^) 



(57a) 



(57b) 



where 



/(to ) = m ( In ^ - 1 



M 



The neutralino mass matrix is 

/ Ml 

Ma 

— M^ cos /3 sin 0^ M^ cos/?cos0^y 

y Mz sin (3 sinOy^^ —Mz sin f3 cos 9 



N 



— Mz cos /3 sin 0^ Mz sin sin 6*^^ \ 

Mz cos P cos 9yy —Mz sin Pcos6^r 

/i 

/X / 



This mass matrix is symmetric and can be diagonalized by a single matrix Z as| 

Z*MnZ~'^ , 



(58) 



(59) 



(60) 



We choose Z to be a real matrix; then the diagonalized neutrino mass matrix can have negative 
entries. We let the entries be CiM^o where M^o are positive masses and takes on a value of +1 or 



-1. The diagonalization can be done numerically, or one can use the analytic expressions 59 

-,1/2 



^1 



e2M^ 



£3 Mo = 



_(ia_ 1^2)1/2 
4 6 ^ 



-{]:a-\c,Yl^ 



(8a-|C2)i/2 
(8a-|C2)i/2 
(8a-|C2)i/2 
3 ' (8a- 1^2)1/2 



-la-lc9- 



2" 



-,1/2 



1/2 



-,1/2 



i(Afi + M2) 



-(M1 + M2) 



-(M1 + M2) 



-(M1 + M2) 



(61a) 



(61b) 



(61c) 



(61d) 
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where 



C2 
C4 



{M1M2 - Mi ~ ^i^) - -(Ml + AhY , 



1 



1 



--{Ml + M2f + -{Ml + M2){MiM2 - Ml - /i") + {Mi + Ma)^^' 
o z 

+ (Mi cos^ 9w + M2 sin^ Ow)MI + ^iMl sin2/3 , 
-(Ml cos^ 9w + M2 sin^ 9w)Ml^i sin 2/3 - MiMa^^ 

+ i(Mi + M2)[(Afi + M2)Ai^ + (Ml cos^ 9^ + M2 sin^ 6'vi/)M| + ^M| sin2/3] 

+ l(MiM2 - M| - /i2)(Mi + M2f - -|;(Mi + M2f , 
16 256 



1 



rRe 



-5 + i(Z3/27)i/2 
1 



21/3 

_4C/3 _ 27S'2 , U 



1/3 



-Co^ - 4C4 



2 ^ 3 ^ 

-C's - — C2 + -C2C4 



(62a) 
(62b) 

(62c) 
(62d) 
(62e) 



These masses given by the above expression are not necessarily such that M^o < M^o < M^o < M^o , 
but the eigenstates can be relabeled. We have corrected a typographical error in the definition of U 
given in Ref. |59(|. The contribution to the minimization conditions is 



AT, 



ix°) _ 



3 r 



AT. 



9^x« 

2 47r2 

1=1 

2^ 47r2 

1=1 



O^^cos/3 + 5^^ sin/3 
g^',sin/3- S*,'^ cos/3 



In- 



In 



M\ 



1 



The factors Q-'j and S*,-^ are defined as||60[ 

Qi'i — [^i3(-^j2 — -^ii tan6'^)] , 
= [Zi4:{Zi2 — ZiitaTi9u])] ei , 



(63a) 
(63b) 



(64a) 
(64b) 



where Ci is the sign of the «th eigenvalue of the neutralino mass matrix. The mixing matrix Z can 
also be given by analytic expressions |p3| 



Zi2 

Zii 
Zii 

ZiA 



1 Ml - e.M^o 
tan6'y[, M2 — eiM^o 



(65a) 



fi[M2 - eiM^o][Mi - ejM^o] - M| sin/3cos/3[(Mi - M2) cos^ By^ + M2 - e^M^o] 



M2[M2 — CiM^o] sin0^;j,[— ^cos/3 + CiM^a sin/3] 



X. 



(65b) 



eiM^o[M2 - eiM^o][Mi - e.M^o] - M| cos^ /3[(Mi - M2) cos^ 6*^ + M2 - CiM^o] 



XiJ 



Mz[M2 - e^M^o] sin6'^[-^cos/3 + e^M^o sin/3] 



and 



Zii 



Zi2 

Zii 



Ziz 
Zii 



Zii 
Zii 



-1/2 



In terms of the mixing matrix Z the bino and gaugino purities are defined as 

Bp = , 



Gi 



Zii + ^12 



(65c) 



(66) 



(67a) 
(67b) 
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respectively. 

The chargino mass matrix is 



Mc - 



M2 \/2Mws,vcl(3 
V2Mw cos f3 -/i 

This mass matrix is not symmetric and must be diagonahzed by two matrices U and 1^ as[j 



where 



U*McV-^ 



\ 0-30+ , detX < ' 



cos (j)± sin 0± 
— sin(/>± cos(j)± 



Here CT3 is the PauU matrix, and 

tan20_ = 2V2A% -J^^^f + ^^^0^ P 



tan2(/)+ = 2V2Mw-r7T o ..,.9 ^ 

^ Af| - /i2 + cos 2/3 



-ficosp + M2 sin/3 



The chargino masses are 



M| + + 2M^ 



± [(M2^ - ^2)2 _^ 4^/4^ cos^ 2P + AM^{Ml + fi^ - 2 Af2/i sin 2/3)] 



1/2 



The contribution to the minimization conditions is 



2 5/Vf^ 



at; 



E 



47r2 



AT, 



2 ffMfi 



47r2 



(^22 COS Sii sin/3 
Qii sin /3 + Sii cos /3 



In — -TT 1 



In- 



- 1 



The factors Qa and Sa are defined as 



Qii — \l ■^ViiUi2 , 



Sii = \/ 2^j2t^ii • 



(68) 

(69) 
(70) 

(71) 
(72) 



(73) 

(74a) 
(74b) 

(75a) 
(75b) 



The Higgs bosons and Goldstone bosons contribute the following contributions in the Landau 
gauge 



gM 



327r2 



2Mw 



cos 6,„ J y" (52 



COS 2/3 ( In _ 1 

M? 
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647r2 cos 6*,^ \ Q"^ 



(76a) 



647r^ cos \ Q 

gMzMl 

647r^ cos 



(- sin 2a + sin 2/?) (in - 1 j . (76b) 



The angle factor a can be eliminated in the above equations using the tree level relations for the 
Higgs masses: 

- 2 cos 2a + cos 2/3 = cos 2/3 ^^^^/^'^/^^^ ^ , (77a) 

2 cos 2a + cos 2/3 = cos 2/3 ( jv^2 _\/-2 ) ' (77b) 

( 2M2 \ 

sin 2a + sin 2/3 = sin 2/3 — ^ , (77c) 

( 2A'P \ 

- sin 2a + sin 2/3 = sin 2/3 — ^ — . (77d) 
The gauge boson contribution is 



35MI 



Then the minimization conditions at one-loop are the following: 

Ti+^AT^^'^ = 0, (79a) 

i 

T2 + ^Ari'' = 0, (79b) 

i 

where i = q,l,lq,ll,x°,X^,H,GB. 
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